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Abstract 

Two central biophysical laws describe sensory responses to input signals. One is a logarithmic relationship between input 
and output, and the other is a power law relationship. These laws are sometimes called the Weber-Fechner law and the 
Stevens power law, respectively. The two laws are found in a wide variety of human sensory systems including hearing, 
vision, taste, and weight perception; they also occur in the responses of cells to stimuli. However the mechanistic origin of 
these laws is not fully understood. To address this, we consider a class of biological circuits exhibiting a property called fold- 
change detection (FCD). In these circuits the response dynamics depend only on the relative change in input signal and not 
its absolute level, a property which applies to many physiological and cellular sensory systems. We show analytically that by 
changing a single parameter in the FCD circuits, both logarithmic and power-law relationships emerge; these laws are 
modified versions of the Weber-Fechner and Stevens laws. The parameter that determines which law is found is the 
steepness (effective Hill coefficient) of the effect of the internal variable on the output. This finding applies to major circuit 
architectures found in biological systems, including the incoherent feed-forward loop and nonlinear integral feedback 
loops. Therefore, if one measures the response to different fold changes in input signal and observes a logarithmic or power 
law, the present theory can be used to rule out certain FCD mechanisms, and to predict their cooperativity parameter. We 
demonstrate this approach using data from eukaryotic chemotaxis signaling. 
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Introduction 

Biological sensory systems have been quantitatively studied for 
over 150 years. In many sensory systems, the response to a step 
increase in signal rises, reaches a peak response, and then falls, 
adapting back to a baseline level, Zsi (Fig. 1 a upper panel). Consider 
a step increase in input signal from /q to I\, such that the relative 

change is Rj = — = ^. There are two commonly observed 

forms for the input-output relationship in sensory systems: 
logarithmic and power law. In the logarithmic case, the relative 

peak response of the system = — = — is proportional 

Z Zst 

not to the input level but to its logarithm Rz ~ log Ri. A logarithmic 
scale of z versus /, namely Az = log /, is often called the Weber- 
Fechner law [1], and is related but distinct from the present 

. . Az , A/ 

definition — = log — . In the case of a power-law relationship, the 
z I 

maximal response is proportional to a power of the input Rz ~ Rj 
(Fig. la lower panel) [2] . In physiology this is known as the Stevens 
power law; the power law exponent /? varies between sensory 
systems, and ranges between /?«0.3 — 3.5 [2]. For example the 
human perception of brightness, apparent length and electrical 
shock display exponents /? = 0.33, 1.1, 3.5 respectively. 

Both logarithmic and power-law descriptions are empirical; 
when valid, they are typically found to be quite accurate over a 



range of a few decades of input signal. For example, both laws 
emerge in visual threshold estimation experiments [3]. In that 
study, the logarithmic law was found to describe the response to 
strong signals and the power-law to weak ones. However the 
mechanistic origins of these laws, and the mechanistic parameters 
that lead to one law or the other, are currently unclear. 
Theoretical studies have suggested that these laws can be derived 
from optimization criteria for information processing [4,5] , such as 
accounting for scale invariance of input signals [6] . Both laws can 
be found in models that describe sensory systems as excitable 
media [7] . Other studies attempt to relate these laws to properties 
of specific neuronal circuits [8,9]. Here we seek a simple and 
general model of sensory systems which can clarify which 
mechanistic parameters might explain the origin of the two laws 
in sensory systems. 

To address the input-output dependence of biological sensory 
systems, we use a recently proposed class of circuit models that 
show a property known as fold-change detection [10,11]. Fold 
change detection (FCD) means that, for a wide range of input 
signals, the output depends only on the relative changes in input; 
identical relative changes in input result in identical output 
dynamics, including response amplitude and timing (Fig. lb). 
Thus, a step in input from level 1 to level 2 yields exactly the same 
temporal output curve as a step from 2 to 4, because both steps 
show a 2-fold change. FCD has been shown to occur in bacterial 
chemotaxis, first theoretically [10,11] and then by means of 
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Author Summary 

One of the first measurements an experimentalist makes to 
understand a sensory system is to explore the relation 
between input signal and the systems response amplitude. 
Here, we show using mathematical models that this 
measurement can give important clues about the possible 
mechanism of sensing. We use models that incorporate 
the nearly-universal features of sensory systems, including 
hearing and vision, and the sensing pathways of individual 
cells. These nearly-universal features include exact adap- 
tation-the ability to ignore prolonged input stimuli and 
return to basal activity, and fold-change detection- 
response to relative changes in input, not absolute 
changes. Together with information on the input-output 
relationship-e.g. is it a logarithmic or a power law 
relationship-we show that these conditions provide 
enough constraints to allow the researcher to reject 
certain circuit designs; it also predicts, if one assumes a 
given design, one of its key parameters. This study can 
thus help unify our understanding of sensory systems, and 
help pinpoint the possible biological circuits based on 
physiological measurements. 



and the power-law exponent /J, depend primarily on a single 
parameter: the steepness (effective HiU coefficient) of the effect of 
the internal variable on the output. 

Results 

Analytical solution for the dynamics of the 11-FFL circuit 
in its FCD regime 

We begin with a common gene regulation circuit [19] that can 
show FCD, the incoherent type 1 feed-forward loop (Il-FFL) [10]. 
In transcription networks, this circuit is made of an activator that 
regulates a gene and also the repressor of that gene. More 
generally, we can consider an input X that activates the output Z, 
and also activates an internal variable Y that represses Z (Fig. 2). 
We study a model (Eq. 1, 2) for the Il-FFL with AND logic (that is, 
where X and Y act multiplicatively to regulate Z), which includes 
ordinary differential equations for the dynamics of the internal 
variable Y and the output Z [20-22]. We use standard 
biochemical functions to describe this system [23]. 

Y = p/{X)-ayY (1) 



dynamical experiments [12,13]. FCD is thought to also occur in 
human sensory systems including vision and hearing [11], as well 
as in cellular sensory pathways [14—17]. 

FCD can be implemented by commonly occurring gene 
regulation circuits, such as the network motif known as the 
incoherent feed-forward loop (Il-FFL) [10], as well as certain types 
of nonlinear integral feedback loops (NLIFL) [11]. Recently, the 
response of an FCD circuit to multiple simultaneous inputs was 
theoretically studied [18]. Mechanistically, FCD is based on an 
internal variable that stores information about the past signals, and 
normalizes the output signal accordingly. We find here, using 
analytical solutions, that simple fold-change detection circuits can 
show either logarithmic or power law behavior. The type of law. 



^ = /4gW^^-«--Z (2) 

The production rate of Y is governed by the input X according to 
a general input function f(X) (in cases where X is a transcription 
factor, X denotes the active state). The maximal production rate of 
Y is Py. The repressor Y is removed (dilution-Hdegradation) at rate 
ay (Eq. 1). We assume here that saturating signal of Y is present, so 
that all of Y is in its active form. The product Z which is repressed 
by Y and activated by X is produced at a rate that is a function of 
both X and Y, denoted G(X , Y). An experimental survey of £. coli 
input functions suggested that many are well described by 
separation of variables: the two-dimensional input function 




input, Rl 



Figure 1. Input/output relationships of sensory systems can be described by a logarithmic law or a power law. a) In many sensory 
systems the dynamical response to a step increase in input signal, /, is a transient increase of output Z followed by adaptation to a lower steady state. The 
relative maximal response is R- = — Two laws are often found. The first is a logarithmic law, R- ~ log i?/. The second law is a power law, R~~R/^' 
with different exponent /J for each sysf'em. b) Fold change detection (FCD) describes a system whose response depends only on the relative change in 
input signal and not the absolute level. Therefore, for a step increase from 1 to 2 and then from 2 to 4 the system response curve is exactly the same. 
doi:1 0.1 371 /journal.pcbi.1 003781. gOOl 



PLOS Computational Biology | www.ploscompbiol.org 



2 



August 2014 | Volume 10 | Issue 8 | e1003781 



Biophysical Laws and Fold Change Detection 



i 

y 

1 




Figure 2. A model for the incoherent feed-forward loop 
includes three dimensionless parameters. In the incoherent type 
1 feed-forward loop (11-FFL) input X regulates an internal variable Y and 
both X and Y regulate Z. The repression of Z by Y is described by a Hill 
function with steepness n and halfway repression point y. 
doi:1 0.1 371 /journal.pcbi.1 003781. g002 

separates to a product of one dimensional functions, 
G(X, Y) = g(X)h( Y) [24] , where g(X) and h( Y) are Hill functions 
(for more explanation see the Methods section). We therefore use a 
general form for the X dependence, g{X), and multiply it by a 
repressive HiU function of Y (Eq. 2), with a maximal production 
rate /J.. The removal rate of Z is a-. Here we consider step input 
functions in which X changes rapidly from one value to another. 
The values of f(X) and g(X) is determined by the step size in 
input. 

For clarity, upper case letters relate to the elements in the circuit 
and lower case letters describe normalized model variables. The 
two-equation model (Eq. 1, 2) has 6 parameters. Dimensional 
analysis (fuUy described in Methods) reduces this to three 
dimensionless parameters (Eq. 3, 4). 

The first parameter, y, is the normalized halfway repression 

ICy- Q 

point of the output, defined by y = where F„ is the pre-step 

Steady state level and Ky- is the level of Y needed to half-way 
repress Z. The second parameter is the cooperativity or steepness 
of the input function, n. The final parameter is the ratio of decay 
a_- „. Y 

0 



— . The normalized variables, y = 



and 



rates of Z and Y, p 
Z 

z = — , are the new dimensionless variables in the 

iS_./a_-^(jr,;<0)' 

model. Table 1 summarizes the parameters in the model for the 
II -FEE. 

This model for the II -EEL describes the response to a step 
increase in input, starting from fully adapted conditions. We 
consider a change between an input level of /o = fiXX < 0) = 1 , to 
a new level I\ =f(X,t>Q) = F. The step is thus characterized by 
the fold change F equal to the ratio between the initial and final 
input levels, F = I\/Iq. 

In order for ECD to hold, the production rate of Z must be 
proportional to X" (g(X)aiX"), where the power law exponent n 
is the same as the HiU coefficient that describes the steepness of the 
input function. In this way, the internal variable, Y, can precisely 
normalize out the fold change in input (see Methods). The model 
thus reads: 



y = F-y 



1 y" 

P 7"+y" 



(3) 



(4) 



The higher y, the more Y is needed to repress Z. The parameter n 
- the Hill coefficient of the input function - is important for this 
study, and determines the steepness of the regulation of the output 
Z by the internal variable Y (Eig. 2). The higher n the more steep 
the repression of Z by Y. The limit «— >oo resembles step-like 
regulation. Biochemical systems often have Hill coefficients in the 
range n=\—4 [23]. The ratio between the removal rates, p, 
describes the relative time scale between Y and Z. For p = 1, Y and 
Z have the same removal rates, and for p» 1, the output Z is much 
faster than Y. 

Goentoro et al. [15] showed, using a numerical parameter scan, 
that this circuit can perform ECD provided that threshold of the Z 
repression, Ky^, is small: that is y^O. We therefore further analyze 
the limit of y— >0, meaning strong repression of Z, where the 
equation for the product Z (Eq. 4) becomes: 



1 



y" 



(5) 



In this limit, the system exhibits fold change detection since it 
obeys the sufficient conditions for ECD in Shoval et al (2010) (see 
Methods). We analytically solved the model (Eqs. 3, 5), in the limit 
of small y, for all values of n, with initial conditions corresponding 
to steady state at the previous signal level, 3^(0) = I , z(0) = Zj, = y" 
(in the limit y— >0). The solution (derived in Methods) is a decaying 
exponential multiplied by a term that contains a Beta function 
(Fig. 3a): 



z(0 



--e-^'(l+p[ 



fF-l 



P F- 



F-\ 



-p,l-«)) (6) 



where the Beta function is i?(a,b,x,y)= J m''' ^{X—uf' ^du. The 

a 

dynamics of the output z shows a rise, reaches a peak Zmax, and 
then falls to the pre-signal steady state (Fig. 3a). At ?-»0 the 
solution is approximately linear with a slope that depends on F, n 
and p: 



■A+p(F"-\)t 



At /— »oo the solution decays exponentially: 



z(t),_ 



■A+e-'" + n 



F-l 



(7) 



(8) 



As in all ECD systems, exact adaptation is found. The error of 

. A z„(F)-zAF=l) „ , 

exact adaptation, £ = — — — goes as e ~ y and 

Zst\F= I) 

vanishes at y— >0. 

We explored how three main dynamical features depend on the 
input fold change F and the dimensionless parameters n and p. 
The first feature is the amplitude of the response, defined as the 
maximal point in the output z dynamics, Zmax. The second 
dynamical feature is the timing of the peak, tpeak- The third feature 
is the adaptation time, T [25,26] which we define as the time it 
takes z to reach halfway between z^ax and its steady state (Eig. 3a). 
We denote Rj as the relative change in the input signal, Rj = F—l 

and R: = -'^ — as the relative maximal amplitude of the 

response. Since p has only mild effects, we discuss it in the last 
section, and begin with p = 1 , namely equal timescales for the two 
model variables. 
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Table 1. A parameter table for the ll-FFL model. 





Parameter 


Biological meaning 


Definition 


ft 
Py 


Mdximdl production r3tG of Y 




«y 


Removal rate of Y 




P. 


Maximal production rate of Z 






Removal rate of Z 




Kyi 


Halfway repression point of Z by Y 




n 


Steepness of input function 






Pre-slgnal steady state of Y 


^/(Z,t<0) 
a.y 


y 


Normalized halfway repression point of Z by Y (dimensionless) 


^ St 


p 


Removal rates ratio (dimensionless) 


a- 



doi:l 0.1 371 /journal.pcbi.l 003781 .tOOl 



A power law relation emerges when the cooperativity n 
of the input function is larger than one; Logarithmic 
behavior occurs when n equals one 

We tested the effects of cooperativity in tlie input function, n, on 
the dynamics of tlie response. Cooperativity seems to have a weak 
effect on the timescales of the response: The adaptation time T and 
the peak time tpeak decrease mildly with the fold F. For p= 1, the 
analytical solution of the time of the peak tpgak for all values of n is: 

tpeuk = log ( , ^' I 1 H ^ I I (see derivation in 

Methods). Substituting the corresponding relative response, R^, 
we receive a mildly decreasing function (Fig. 3b). 

In contrast to the mild effect of cooperativity on timescales, 
cooperativity has a dramatic effect on the response amplitude. The 
maximal amplitude of the output z relative to its basal level, R:, 
increases with the fold and behaves differently for each n. For low 
steepness, n=l, R~ increases in an approximately logarithmic 
manner with Rj (for R; » 1), R. ~ \og(F — 1) = log R/ (normalized 
root-mean-square deviation, RMSD = 0.0\4 for fitting to 
alogRi + h compared to RMSD = 0.029 for fitting to aR/'- see 
Methods). More precisely the analytical solution is 
R-=product\og^ (see Methods) (Fig. 4a). The function 
y= product \og{x) is defined as the solution to the equation 
x=ye^. The producdog function is approximately linear at x«l, 
and approximately log(.x:) at x » 1. 

For n = 2, the peak response increases linearly R^^F —l=Rf. 
For w = 3, the increase is approximately quadratic, 
Rz~{F-\)^ = Rl^ (Fig. 4b). We find that for any «>1, the 
increase is approximately a power law with exponent n — l in the 

limit of large Rf. R,c^ -^(F -l)"'^ (see Methods) 
(i^MSi) = 0.049,0.003,0.009 for fitting to ^R/ compared to 

RMSD = 0.m, 0.436, 0.671 for fitting to a\ogRi + b for 
n = 2, 3, 4 respectively). Note that the pre-factor in the power 
law is also predicted to depend simply on the HiU coefficient for 

«> 1, namely to be equal to j- (for p= 1). Indeed in fitting the 

numerical solution the best fit parameter is approximately 
P~n-\: jS = 0.92, 1.99, 3.0 for n = 2, 3, 4 respectively. The 
dependence of output amplitude on input fold-change is thus a 



power law, similar to Stevens power law, except for n = 1 where 
the output dependence is logarithmic. 

One point to consider regarding step input functions is that 
realistic inputs are not infinitely fast steps; however, a gradual 
change in input behaves almost exactly like an infinitely rapid step, 
as long as the timescale of the change in input is fast compared to 
the timescale of the Y and Z components. To demonstrate this, we 
computed the response to changes in input that have a timescale 
parameter 0£.v that can be tuned to go from very slow to very fast: 
jr(0=l + (i^-l)tanh(o(.^/) (Fig. 5a). When a^> ~o(,,,o(-, the 
behavior of the relative maximal amplitude of the response, R^, 
as a function of the relative change in the input signal, Ri = F —I, 
is very similar to the infinitely fast step solution (less than 5% 

ax 

difference for F = 50, n = p= \ and — = 10, Fig. 5b). When the 

change in input is much slower than the typical timescales of the 
circuit, the response is very small, since the signal is perceived 
almost as a steady-state constant. For slow changes in input, the 
Il-FFL response can be shown to be approximately proportional 
to the logarithmic temporal derivative of the signal [27-30]. 

A nonlinear integral feedback mechanism for FCD also 
shows a power law behavior 

In addition to the Il-FFL mechanism, a non-linear integral 
feedback based mechanism (NLIFL) for FCD at small values of y 
has been proposed by Shoval et al [1 1] (see Methods section) 
(Fig. 6a). This mechanism is found in models for bacterial 
chemotaxis [28]. The fuU model is described by: 

Y = kY(Z-Zo) (9) 

^ = i9.-gm^^;^-^zZ (10) 

Its dimensionless equations following dimensional analysis (fully 
described in Methods) are: 

i'=X-^) (11) 
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tpeak 

Time 



(b) 



Time of the peal< 



Adaptation time 



05 
<D 
Q. 





Figure 3. The 11-FFL shows FCD in the limit y->0. a) Response to a step increase in input from /o to Ii, which can be described by 

Rj s — - =F— 1 where F= ^. The output dynamics show three features of interest: the amplitude of the peak response r^ax the timing of the 

peak, tpi,,ik and the adaptation time t. b) The time scales of the response, the timing of the peak tp^ak and the adaptation time t, mildly decrease with 
the relative change in the input signal, Rj. The steepness, n, does not have a dramatic affect on this decrease. 
doi:1 0.1 371 /journal.pcbi.1 003781. g003 



--F" 



r 



1 . 



Where the new variables are: y 



dimensioiiless parameters are 



Y 

defined 



Z 

as: ■ 



(12) 



and the 



and 



■ (Methods). Table 2 summarizes the parameters 



in the model for the NLIFL. 

We solved the NLIFL model (Eqs. 11, 12) numerically for the 
limit 7« 1 and find that the maximal response increases with the 
relative change in the signal in a power-law manner, R- = al^\ 
(Fig. 6b). The best-fit power law exponents increase with «, 
namely iS = 0.7, 1, 1.4, 1.8 at M=l,2, 3,4 for p=\. A \ogRi 
dependence does not fit the data at n 
to aRi^^ compared to i?M5'Z) = 0.1 8,0.39,0.55,0 



1 [RMSD = ^a\ for fitting 
for fitting to 



fllogi^/ for «= 1, 2, 3, 4 respectively). To a good approximation, 
the power law is linearly related to the steepness parameter w, by 
j5 = 0.1+0.46« (Fig. 6c). 

The time scales in this circuit seem to decrease faster with the 
fold F for n>2 than in the Il-FFL case, tpeak,T: — cR'} where 



-0.3,-0.6,-0.9,-1.3 and e,-- 



1,-1.5 



-0.3,-0.7, 

at «= 1, 2, 3, 4 (Fig. 6d, all the fits of tpeak,T: have R^>0.99). 

Given the results so far, one can use the present approach to 
rule out certain mechanisms. If one observes a logarithmic 
dependence, one can draw at least two conclusions: (i) the NLIFL 
model addressed here can be rejected, (ii) if the Il-FFL model 
addressed here is at play, its steepness coefficient is w = 1 . 

If one observes a hnear dependence of input on output, the Il- 
FFL and NLIFL mechanisms cannot be distinguished. The 
steepness can be inferred to be about n = 2 for both circuits. 

Logarithmic law in eukaryotic signaling FCD 

We applied the present approach to data from Takeda et al 
[17] on Dictyostelium discoideum chemotaxis. In these exper- 



PLOS Computational Biology | www.ploscompbiol.org 



5 



August 2014 | Volume 10 | Issue 8 | e1003781 



Biophysical Laws and Fold Change Detection 




(b) 



N 2 




100 



200 



300 



400 



500 



Ri 



Figure 4. The response amplitude follows an approximately logarithmic law for /7= 1 and a power law at n>^ . a) The numerical solution 

of the amplitude of the response, R-, with n = ^,\s shown (blue dots) as well as its analytical solution R-= product log — (blue curve) in a log-linear plot. 

e 

A fit to R: = alogRi+b (red curve) captures the behavior better than a fit to R: = aRi'' (green curve), b) For n>1 the numerical solution of the 

amplitude of the response, R^, is shown (in dots) as function of Rj in a log-log plot. A fit to a power law behavior R^ = ^^/'' °"'y parameter 

(solid lines) describes the numerical results better than a fit to a logarithmic behavior R:=alogRi+b (dashed lines). At the limit of large Rj, jl~n — l. 
doi:1 0.1 371 /journal.pcbi.1 003781. g004 



iments, the input is cAMP steps applied to cells within a micro- 
fluidic system, and the output is a fluorescent reporter for Ras- 
GTP kinetics. The output showed nearly perfect adaptation 
and FCD-like response to a wide range of input cAMP steps. 
We re-drew the peak amplitude (Fig. 7a) and the time of peak 
(Fig. 7b) as a function of the added cAMP concentrations and 
find that it is well described by the analytical solution of the 



maximal response and time of peak for an Il-FFL circuit with 
n=\. The peak amplitude [R^] as a function of the relative 
input R] is well described by a logarithmic relationship (mean- 
square weighted deviation, MSWD = Q.Q9 for fitting the data 
to a\ogRi-\-b considering the error-bars - see Methods). 
Fitting it to a power law uR/^ results in a small exponent 
jS~0.12 (MSWD = Q.\\) (Fig. 7c). Such a small power law 
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(a) The relative response for n=1 



N 

a: 




5 10 50 100 



Ri 

(b) The relative response for n=2, 3, 4 




5 7 10 15 20 30 

Ri 



Figure 5. Rapidly cKianging input signal leads to responses similar to a step increase in signal; slowly changing input leads to weak 
response, a) Input signal with a tunable timescale, X(<)= 1 + (^— l)tanh(5(vO with F=10, = This signal goes form level 1 to level f, with a 
halfway time that goes as Xja.^. b) The relative maximal amplitude, R-, as a function of the relative change in the input signal R[, is plotted for various 
values of the input timescale a,. When the signal changes much faster than the timescale of the circuit, s(v»o;|.,o!_-, the response is similar to the 
analytical solution for an infinitely fast step in input (Full red curve). When the timescale is slow, a.x<^n.y,^-, the response of the circuit is weak. 
doi:1 0.1 371 /journal.pcbi.1 003781. g005 



exponent can only be obtained with a negative cooperativity in 
the NLIFL model considered here. Such negative cooperativity 
is rare in biological systems [31,32]. If we consider only 
positive cooperativity («>1), as found in most biological 
systems, the NLIFL model considered here provides a poor 
fit to the data (j5 = 0.7,M5I^'Z) = 0.17) (Fig. 7c). 

Thus, the present analysis is most consistent with an Il-FFL 
mechanism considered here with n=\. The same is found when 
plotting the observed time-to-peak it peak) versus the analytical solution 



of the Il-FFL model (?, 



'peak = log ( ( 1 + ^ 

(i + i^_-)'"-i 



witli 



[MSWD-- 



Rr 



l+Rl 



\ 



■0.02 



for 



fitting 



to 



a log 



V 



product log — , 



(Fig. 7d). This agrees with the 



numerical model fitting performed by Takeda et al, who conclude 
that an 1 1 -FFL mechanism is likely to be at play (they used « = 1 in 
their 1 1 -FFL model, which is based on degradation of component Z 
by Y, rather than inhibition of production of Z by Y as in the present 
model). 

In this analysis we assumed that the experimentally measured 
fluorescent reporter is in linear relation to the biological sensory 
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(a) 



(b) 



The relative response for n=1 , 2, 3, 4 




5.0 7.0 10.0 15.0 20.0 
Rl 

(<=) (d) Tinnesclalesforn=1,2, 3, 4 




The steepness, n 5 10 15 ' '° 

Rl Rl 



Figure 6. A different circuit showing FCD, tKie non-linear integral feedback loop (NLIFL), also exhibits a power law behavior, a) The 

NLIFL mechanism, b) The amplitude of the response is a power law of the relative change in input signal, c) The power-law exponent (i increases 
linearly with n. d) The time-scales decrease faster with the fold change of the signal, Rj, and with n than in the incoherent feed-forward loop case 
(Fig. 3b). 

doi:1 0.1 371 /journal.pcbi.1 003781 .g006 



output, Ras activity. If this relation turns out to be nonlinear, the 
conclusions of this analysis must be accordingly modified. 

Effect of timescale separation between the variables 

In the eukaryotic chemotaxis system, the two model variables Y and 
Z have similar timescales (p = '^-/^^ ~ 1). We also studied the effect of 
different timescales (p # I), and find qualitatively similar results. A 



logarithmic dependence of amplitude on F is found when n=l, and a 
power law when n>l. The power law p increases weakly with p 
(Fig. 8a). In the limit of very fast Z (p^ co), the solution approaches an 
instantaneous approximation (obtained by setting z = 0) in which the 
power law is n instead of n — 1 (Fig. 8b). There is a cross over from the 
Stevens power law /? = w — 1 when F»p, to the instantaneous model 
power law j8 = « when 1 « _F « p (Fig. 8c). An analytical solution that 
exemplifies this crossover can be obtained at «=1, where 



Table 2. A parameter table for the NLIFL nnodel. 





Parameter 


Biological meaning 


Definition 




k 


Zo 


Steady state level of Z 








Maximal production rate of Z 








Removal rate of Z 






Ky 


Halfway repression point of Z by Y 






n 


Steepness of input function 






yO 


Pre-signal steady state of Y 


^ //lg(Z,t<0) ^ 


t 


y 


Normalized halfway repression point of Z by Y (dimensionless) 






p 


Timescale ratio (dimensionless) 


klS-^g(X,t<Q) 





doi:l 0.1 371 /journal.pcbi.1 003781 .t002 
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R- = pproduct \og( — ) (Methods). Because of the limit behavior of 
pe 

the productlog function mentioned above, at small fold values R^^Rj, 
and at large values R^^logRj. In summary, the instantaneous 
approximation, commonly used in biological modeling, must be done 
with care in the case of FCD systems. 

Discussion 

This study explored how two common biophysical laws, 
logarithmic and power-law, can stem from mechanistic models of 
sensing. We consider two of the best studied fold-change detection 
mechanisms, and find that a single model parameter controls which 
law is found: the steepness n of the effect of the internal variable on 
the output. We solved the dynamics analytically for the Il-FFL 
mechanism, fmding that logarithmic-like input-output relations 
occurs when n=l, and power-law occurs when n>\, with power 

law P = n—\, and prefactor j- at p = 1 . The nonlinear integral 

feedback loop (NLIFL) mechanism - a second class of mechanisms 
to achieve FCD - can only produce a power law. Thus, if one 
observes logarithmic behavior, one can rule out the specific NLIFL 
mechanism considered here. This appears to be the case in 
experimental data on eukaryotic chemotaxis [17], in which good 



agreement is found to the present results in the 1 1 -FFL mechanism 
with n = 1 in both peak response and timing. 

This theory gives a prediction about the internal mechanism for 
sensory systems based on the observed laws that connect input and 
output signals. Thus, by measuring the system response to different 
folds in the input signal one may infer the cooperativity of the input 
function and potentially rule out certain classes of mechanism. For 
example, if a linear dependence of amplitude on fold change is 
observed (power law with exponent /J= 1), one can infer that the 
steepness coefiicient is about n = 2 for botii the specific Il-FFL and 
NT TFT , circuits considered here, with slight modification if the 
timescales of variables are unequal. Such a linear detection of fold 
changes may occur in drosophila development of the wing imaginal 
dLsk [33-35]. 

The problem of finding the FCD response amplitude shows a 
feature of technical interest for modeling biological circuits. In many 
modeling studies, a quasi-steady-state approximation, also called an 
instantaneous approximation, is used when a separation of timescales 
exists between processes. In this approximation, one replaces the 
differential equation for the fast variables by an algebraic equation, by 
setting the temporal derivative of the fast variable to zero. This 
approximation results in simpler formulae, and is often very accurate, 
for example in estimating MichaeUs-Menten enzyme steady states [36] . 
However, as noted by Segel et al [36], this approximation is invalid to 
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Figure 7. A mechanism of euicaryotic signaling FCD illustrates this theory, a) The response of Ras-GTP to different concentrations of added 
cAMP in Dictyostelium discoideum chemotaxis is re-plotted together with the timing of the peak (b). A logarithmic function describes the data well. 
The black lines are our fit to the data, c) The response of Ras-GTP is re-plotted as function of the different fold changes in cAIMP concentrations. The 
solid line is a fit to ^ a log Rj + h, the black dashed line is a fit to R^ ^ciRj^' and the green dashed line is a fit to a power law with exponent /i = 0.7. 
d) The corresponding solution for the timing of the peak for Il-FFL with n = 1 explains well the data. 
doi:1 0.1 371/journal.pcbi.1 003781 .g007 
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Figure 8. The instantaneous approximation does not capture the correct amplitude behavior, a) The power law for n = 1 increases mildly 
with p to a value between 1 and 2. b) The instantaneous approximation (in red) and the full model solution (in black) are plotted as function of time 
for i^ = 5,« = 2,p 3:250. c) The maximal response r^ax normalized to ;'™'' = i^" is plotted for different folds and for n = 1. The error between the 
maxima of the instantaneous approximation and the full model increases with the fold F. 
doi:1 0.1 371/journal.pcbi.1 003781 .g008 



describe transients on the fast time scale. In the present study, we are 
interested in the maximal amplitude of the FCD circuits. In some input 
regimes, namely F » p, the instantaneous approximation predicts an 
incorrect power law. To obtain accurate estimates, the full set of 
equations must be solved without setting derivatives to zero. 

It would be interesting to use the present approach to analyze 
experiments on other FCD systems, and to gain mechanistic 
understanding of sensory computations. 

Methods 

The two dimensional input function can be considered 
as a product of one dimensional input functions 

Consider a general partition function for an input function with 

aX" 

an activator and a repressor: G(X, Y) = ; . 



The regime in which FCD applies is that of strong repression, 
yY'"»l and non-saturated activation /JX" « 1 [10]. In this 

limit, G{X, Y) ~ , and is thus well approximated by a product. 

yym 

More generally, G(X,Y) is a product of two functions 

aX" 1 

whenever binding is independent, G(X, Y) = 



l+ZJX" i+yy™ 

, which occurs when the relation 3 = 



aX" 

\+pX" + yY'" + liyX"Y" 
holds. The biological meaning of the relation is that X and Y bind 
the Z promoter independendy so that the probability of X to bind 
the promoter and the probability of Y to unbind equals the 
multiplication of the probabilities: 

P(Xbound A Yunbound) = P(Xbound)P(Yunbound). 
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In the NLIFL case, one can show from the MWC model 
chemotaxis by Yu Berg et al [28] that in the FCD regime it is 
simply a power law. 

Dimensional analysis of the full model for 11-FFL and 
NLIFL 

We performed dimensional analysis of the full model of the II- 
FFL (Eq. 1, 2) by rescaling as many variables as possible. The 
rescaled variables: 



y= 



(Ml) 



l}Ja,g{X,t<0) 

Where Y° is the pre-signal steady state of Y, derived by taking 

y = 0: F2=^/(Jr,t<0), and Z,,= ^g(jr) ^ is the 

Steady state of Z derived by taking Z = 0. Substituting these 
rescaled variables we receive: 



J^^ f{X) Y_ 

yo '■yf(x,t<Q) '"yo 



(M2) 



Z g{X) K"y^ Z 

fSz/oizg{X,t< 0) g{X,t< 0) k;^ + Y" Pz/o^^g(x,t < 0) 

Since we assume that /{X) is determined by the step size in input, 

f(^X) 

we can consider merely the fold change F in input, F= — — — . 

For FCD to hold we consider F" = -—^ — ^— . Defining the 

K-yz 

rescaled repression threshold y= — ^ we receive in the new 
rescaled variables (lower case letters y and z): 

—y = F-y (M3) 



1 y" 

JL2 = F"— z 



Rescaling the time to t-*oiyt and defining p = — yields to Eq. 3, 4 

^y 

in the main text. 

We also performed dimensional analysis of the fuU model of the 
NLIFL (Eqs. 9, 10) by rescaling as many variables as possible. The 
rescaled variables: 



Where Yf, is the pre-signal steady state of Y, derived by taking 

7 n A ■ 7 7 yO V fMlA^ lV'" A 

Z = 0 and assuming Zy, = Zq : i^, = Ky I — 1 ) , and 



Ky /Pzg(X,t<0) 



^ yo V o!-^o 

iables we receive: 



-1 



KzZn 



Substituting these rescaled var- 



1 f Y _ Y Z f 



kZo 1 +7" yo YO,^Zol+y" l+r 



Z f _,z f ^ 



(M5) 



Zol+>'« Zol+y" K^,z+Y" -l+y"Zo 

After algebraic manipulation and in the new rescaled variables 
(lower case letters y and z): 



1 y" 
kZol+r 



y=y{z- 



1-1- y"^ 



(M6) 



1 gW / 

az g(X,r<0)y«+>'" 

gW 



We consider here also ^ 

g(A',/<0) 

Rescalmg the time to t^kZQ — - — 1= 



defining 



kpzgiyj<o) 



kt and 

y" tXz 
yields to Eqs. 11, 12 in the main text. 



Proof that FCD holds in the model for 11-FFL and NLIFL 

Given a set of ordinary differential equations with internal 
variable y, input F and output z: 



j=/(y,z,F) 



2=g'(y.z.F) 



(M7) 



(M8) 



According to Shoval et. al. (2010), FCD holds if the system is 
stable, shows exact adaptation and g and f satisfy the following 
homogeneity conditions for any p>Q: 



f(py,z,pF)=pf(y,z,F) 



g(py,z,pF)=g(y,z,F) 



(M9) 



(MIO) 



In the model for 11-FFL (Eq. 3, 4) at the limit of strong repression 
y-»0: 



y= 



(M4) 



f{py,z,pF) = (pF -py) =pf{y,z,F) 



Z f 

"zoi+r 



1 

-g{py,z,pF\^Q'^(pFf -z = F" — -z = g(y,z,F\^f^ 
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Exact adaptation also holds at y->0, Zst—zo- This holds also for 
the NLIFL (Eqs. 9, 10). 



the amplitude of the maximal response, /(7?z ,7?/ ,«,p) = 0, yielding 
an intractable equation: 



Analytical solution for the 11-FFL 

The solution for y is an exponent: 

y{t) = e-'+F{\-e-') 



(Mil) 



The general solution for the ODE z=f(i) — z with the initial 
condition z(0) = zo is: 



z(0 = zoe"' + e"'|ey(T)(/T (M12) 

0 

For our model Eq. Ml 2 reads: 

z(0 = /e-'+/'pe-'|e^F"^-3,^^^l— (M13) 



By changing the variable in the integral in Eq. Ml 3: w = e^ 
we get: 



z(0=A-'"(i+p(-^ 



F 



-p-1 



il-w)-"dw) (M14) 



Which is by definition the solution in Eq. 6. 

Analytical solution for the time of peak for the ll-FFL 

At the time of peak z = 0, therefore from Eq. 5 in the main text 
we get: 



'{t peak) J 



(M15) 



From our definition of the relative response = ^ ^""^^ — we 



have: 



y{tpeak) = 



X- + 1)'" 



(M16) 



Substituting the solution of y (Eq. Mil) and by algebraic 
manipulation we receive the analytical solution for tpeak'- 



imax=log 



Ri 



\+R,\ (1+ii^ 



(M17) 



Analytical solution for the maximal response 

The analytical results were derived by taking the derivative of 
the solution for z{t) (Eq. 6 in the main text) and substituting time 
dzitpeak) rv 

of the peak (Eq. Ml 7), — -j =0. This provides an equation for 



f{R,,Ri,n,p) = I ^-^] ' - (1 + i^,) ( 1 - ^ 



Ri 



{l+R.f"^ 

Where we used the identity: 

B{a,b,-pA-n) = i-iy'B{b-\a-\p + n,l-n). This identity 

1 

can be easily proven by using the change of variable, w= in 

u 



the integral of the Beta function. 
For n = 1 Eq. Ml 8 becomes: 



Using the Series function of Mathematica to expand Eq. Ml 9 in 
the limit of large Rj, R^ and keeping high orders in Rj, R^ yields: 



l-(l + ;^)'(l+^z)-plog(^^j=0 (M20) 
Using [iH — j ~1h — in the limit of large x we receive: 



-/?,-p-plog(^|^j=0 
Taking the exponent of this Eq. M2 1 yields: 



pe p 



(M21) 



(M22) 



The solution for Eq. M22 is by definition the productlog function: 
Rz = p product log 1 — I . 
For /9= 1 Eq. M18 becomes: 



Ri ' R. 



b 

Since B{a,b,2S!))= |y^'^-^= " C'' " ^) " ' ^q- M23 



yields: 

l+R 



By algebraic manipulation Eq. M24 becomes 1 + _Rz + 
log ( ^ ) = 0- Taking the exponent of this equation yields: 



Ri 



-R.e'' 



(M25) 
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The solution for Eq. M25 is by definition the productlog function: 
J?z = product log \' 



For n>\,p=\ we define f = (1 — 1, substituting this 
new variable into Eq. M18 we have: 



^-(l + i)"^+(-l)"5(^,^,l+«,l-«)=0(M26) 

Using the Series function of Mathematica for large Ri and c 
yields: 



(«-l)i^/«-^ + (l + O"+') = 0 



(M27) 



Keeping the highest order in Rj and £, we receive: 
in-\)C^Ri"~^. RecaU that f « (1 + {)" = 1 + J?^ « R, for large 

Rz and and therefore i?z— r^/" '- 

« — 1 



The instantaneous approximation does not capture the 
correct amplitude behavior 

For the instantaneous approximation to be true at large p, the 

error, (Fig. 8a), between the maximal amplitude in the 

instantaneous approximation and the full model should vanish at 
p-»oo. 



Fits and numerical simulations 

All the numeric simulations and fits were made in Mathematica 
9.0. 

The root-mean-square deviation (RMSD) [37] calculated for 
comparing the goodness of fit between the two models is defined 

1 



as: RMSD = 



'-^i, numeric -^ijinodel* 



The data points from Takeda et al were extracted 1)y using the 
'ginput' function of MATLAB. The fits for the data were made 
using the NonlinearModclFit function considering the error-bars. 

Ax,-, as weights, w,- = - — ^. 

The goodness of fit was tested using the mean-square weighted 
deviation (MSWD) [37] which sums the residuals (r) - sum of 
squares of errors with weights of \/s.d.: MSWD = 
1 



1 " 
n— \ ■f— f 



Note on biophysical law terminology 

Az 



A/ 



We define logarithmic response as — = log — . In contrast, 

traditional definition of the Weber-Fechner law (also called the 
Fechner law) in biophysics is (e.g. ref [3]) as Az~log/. Thus the 
present definition concerns relative change in input and output, 
whereas the Weber-Fechner law concerns absolute input and 
output. Note also that the Weber-Fechner law is distinct from 
Weber's law, on the just noticeable difference in sensory systems, 
whose relation to FCD was discussed in Ref [1 1]. 



Az,v 



1 



tpeak ~{F- \)tpeak{F) (M28) Acknowledgments 



Where tpeak{F) decrease with F slower than F^^ , therefore 
Az,„5r 

'— ~f(F) with/(i<) a monotonic increasing function of-F. This 

proves that even at large p, the error increases with F (Fig. 8b) and 
can be very large. 
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